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Abstract 

The Fisher-KPP equation is a model for population dynamics that has generated a huge amount of 
interest since its introduction in 1937. The speed with which a population spreads has been computed 
quite precisely when the initial data decays exponentially. More recently, though, the case when the 
initial data decays more slowly has been studied. Building on the results of Hamel and Roques TO, 
in this paper we improve their precision for a broader class of initial data and for a broader class of 
equations. In particular, our approach yields the explicit highest order term in the location of the level 
sets. In addition, we characterize the profile of the inhomogeneous problem for large times. 


1 Introduction 


In this article, we study the long time behavior of solutions to the Fisher-KPP equation 

\ ut = Uxx + t>0,xGM., 

[ it(0, x) = X G K., 


(I.l) 


where / is an L-periodic function in x, for some L > 0, which vanishes when u is 0 and 1 and where uq is a 
function taking values between 0 and I that decays slowly as x tends to infinity. This equation comes from 
adding a diffusive term to the logistic equation, and it was originally introduced, with f{x,u) = u{l — u), 
as model for population dynamics in the first half of the twentieth century [11,18]. In this interpretation, u 
stands for the local population density normalized by the carrying capacity, while / accounts for the effect 
of an inhomogeneous environment on the growth of this population up to the carrying capacity. 

Fisher and Kolmogorov, Petrovskii, and Piskunov showed that if the initial data is localized then the 
population would spread at a constant rate in the form of a traveling wave. This has been made more precise 
and succinct over time [5,6,13-15,17,22,25,26]. In these works, the highest order terms in the spreading 
rate are obtained by studying the linearized problem. When uq is not localized, but decays exponentially, 
like the behavior is similar. Here, the population converges to a traveling wave, the speed of which is 

determined by the order of the exponential decay and is given again by the linearized problem [26]. 

Recently, the case when uo decays slower than any exponential and / is homogeneous was investigated by 
Hamel and Roques [16]. They show that u becomes asymptotically flat in time, which we will state precisely 
later. This implies that a traveling wave cannot give the limiting behavior of u. As a result, they investigate 
the location of the level sets of u and show that, for each m G (0,1), the level sets satisfy 


{x G M : u(t, x) = m} C 




( 1 . 2 ) 


where Cm < Cm are constants depending on m. These results hold for only a small class of initial data; 
for more general data, they show that level set is bounded above and below by (exp —(/'(O) ± e)t) for 
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any e > 0. The main results of this paper are twofold: firstly, in the homogeneous setting, we explicitly 
characterize the constant in (1.2) for a broad class of the admissible initial data; secondly, we obtain equally 
precise results in the previously unstudied inhomogeneous setting. In fact, in the inhomogeneous setting, we 
show that u grows at each point according to the solution to a generalized logistic equation. This provides 
an analogue to the asymptotic flatness result of Hamel and Roques in the inhomogeneous setting. 


2 Main results 

We now state our results more precisely. First we assume that there is a (5 > 0 such that f{x, u) is a positive, 
function with respect to u that satisfies 

f{x, u) < fu{x, 0)u for all u G [0,1], 
is decreasing in u, and 

u B ’ 2.3) 

/„(x, 0) > 0 and 

/„(x,1) < 0 

We note that some of the assumptions on the signs of / and its derivatives can be relaxed, but we use them 
here to simplify the arguments. 

We assume that Uq is a positive function that decays to zero slower than any exponential. By this we 
mean that, for any e G (0, oo) 


lim 'Uo(a:) = 0, and lim 'Uo(a;)e'^“ = oo. (2.4) 

X—>-oo X—>-oo 

For instance, functions which decay algebraically as x tends to infinity satisfy the conditions in (2.4). Though 
our technique can handle any functions of this type, the presentation is much simpler if we restrict to initial 
data, uo, which are eventually monotonic and which decay slower than for any /3. Namely, we assume 

that there is a point xo > 0 such that 


Uq{x) < 0 for all x > xq, and 


lim 

X—>-oo 


Uo(x) 

uo(x)/log(uo(x)) 


= 0 . 


(2.5) 

( 2 . 6 ) 


We will comment on the effect of this assumption later. In addition, since we are studying the effect of a 
slowly decaying tail at oo, we additionally assume that uq is front-like. In other words, we assume that 


liminfMo(a^) > 0. (2.7) 


Again, we note that the positivity assumption on uq can be relaxed, but we use it here to simplify the 
presentation. 

Classical results address the case when uq ~ [5,6,11,18,22,26]. If A > 1 then u converges to a 

traveling wave with speed 2 in the correct moving frame, while if A < I, it converges to a traveling wave with 
speed A -b I/A in the moving frame. In short, this means that the level sets of ito move with a constant speed, 
which increases to infinity as A decreases to zero, and that any two level sets travel at the same speed and 
remain a bounded distance from one another for all time. Hamel and Roques point out in [16] that, using 
these as a family of sub-solutions, we can see immediately that the level sets of any solution with initial data 
like eq. (2.4) will spread to the right super linearly in time. We now calculate the speed of these level sets 
explicitly. 


2.1 The homogeneous case 

Our first result concerns the spreading of the level sets in the homogeneous setting. 
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Theorem 2.1. Suppose that u solves eq. (1.1) with f depending only on u and not on the space variable, x. 
Suppose further that f satisfies condition (2.3) and that uq satisfies conditions (2.4) and (2.7). In addition, 
let Uq satisfy conditions (2.5) and (2.6). Let ip be the unique (up to translation) global in time solution to 


‘Pt = f{<p)- 

Then, for each m € (0,1), define Tm € M. to be the unique time so that p{Tm) = rn. Then, for any r > 0, 
we have 

limsup sup u{T,x)<m. (2.8) 

T^oo x>tt”^(v3(Tm-T))-rT 

and 

liminf inf u(T,x)>m. (2.9) 

T-s-oo x<«”^(v3(Tm-T))+rT 

In order to make this result clearer, we first consider a specific example. If / is the standard Fisher-KPP 
non-linearity, i.e. f(u) = u(l — u), then 

T 

<f{T) = ^ implying that ip{Tm -T)= ^ e~'^ -k 0(e~^'^), 

1 -|- e 1 — TO 

where 0(e“^^) is bounded by a constant multiple of for all T. For a general non-linearity, /, one 

can easily see (by linearization) that, up to leading order, p{Tm — T) = Cme~^ W’"^ where Cm is a constant 
depending only on / and m. In addition, defining the level set of u as 

Emit) {x : uit,x) = to }. 


if we look at, e.g., initial data such that u^ix) = x “for large enough x, theorem 2.1 implies that 


Emit) 



( 2 . 10 ) 


Notice that the non-linearity affects the constant as the linearized equation would yield only 1 /to. The 
fact that the full non-linearity affects the highest order behavior is in contrast to the case when uq decays 
exponentially, where the linearized equation completely determines the highest order behavior. 

An easy corollary of theorem 2.1 is the following. 

Corollary 2.2. For any e > 0, under the assumptions of theorem 2.1, there exists Tq > 0 depending on all 
parameters such that if T > Tq then 

EmiT) C [uq ^ ipiTm+e “ T)),u0^ ipiTm-e “ T))] . 


Again, it is instructive to look at the example that we considered above. Corollary 2.2 implies that, for 
any e and any t large enough, we have 


Emit) C 



Remark 2.3. We comment briefly on the assumption on the decay of the first derivative ofuo, eq. (2.6). In 
practice, this means that uq decays slower than e~^. We could remove this restriction, assuming only that 
Uq satisfies eq. (2.4). In this case, we would obtain a linear in time error term in eqs. (2.8) and (2.9) and in 
corollary 2.2. In other words, there is a constant r > 0 such that, for any to G (0,1) and e > 0, we have 

EmiT) C [uq^ ipiTm+e “ T)) - vT, u(f^ fpiTm-e “ T)) + vT] , 

for T large enough. We omit proving this here, in the interest of presenting a single, clear argument and 
result. 


3 









2.2 The inhomogeneous case 

When the non-linearity, /, depends periodically on x, the level sets fail to capture the spreading behavior. 
Intuitively this is because the function uq is very flat at infinity, so it will grow as if it were a solution to 
the periodic, inhomogeneous Fisher-KPP equation. In particular, u will oscillate in space enough that the 
width of the level sets will grow to infinity as t tends to infinity. In other words, the distance between the 
rightmost and leftmost points of the level set Emit) tends to infinity. 

To circumvent this, we look instead at the locations of the intervals of length L, the length of the period 
of /, of a given mean oscillation. In other words, we seek to describe the location of the average-level 


Emit) 




where we define, for any interval [a, b] and any function 7 ( 0 :), 



Using Emit), we find similar results to theorem 2.1 in the inhomogeneous setting. 

As we alluded to above, instead of the logistic equation, the behavior of u on each interval of length L is 
captured by the global in time solution to 


\'ft ='^^xx+ fix,ip), for all (t, x) € R X T, ^ 

\ffiO,x)dx = i, 

where T = [0, L] with the ends identified. In other words, ip is an L-periodic solution to the equation 
above. We briefly note that, although the techniques to show the well-posedness of eq. (2.11) are similar to 
those used in the theory of well-posedness for the Fisher-KPP equation, see e.g. [2], we are unable to find a 
treatment of it in the literature. As such, we dedicate section 6 to proving the existence and uniqueness of 
the solution to this problem. 

We define one last quantity before we state our results. For each m G (0,1) and each T > 0, we may 
solve the terminal value problem 

r + fix, ipE-^), on [0, T] X [0, L], 

S T fo^ ‘p'^’’^iT,x)dx = m, (2.12) 

[ipE--i0) = Bim,T), 


with periodic boundary conditions and where Bim,T) is an unknown number depending only on /, m, and 
T. It is simple to check that the system eq. (2.12) is well defined and is continuous in m and T. It is also 
easy to check that i3(m,T) is increasing in m and decreasing in T. 

We now formulate our results precisely. 

Theorem 2.4. Suppose that u solves eq. (1.1) and that f is an L-periodic function satisfying the conditions 
in eq. (2.3). Suppose uq satisfies the conditions in eqs. (2.4) to (2.7). For any r > 0, we have that 


rx+L 

limsup sup -r uiT,y)dy<m, 

T->oo x>u~^{B{m,T))-rT Jx 


and 


rx+L 

liminf inf -r uiT,y)dy>m, 

T->-oo x<u-^{B{m,T))+rT JX 


where Bim,T) is defined in (2.12). 


(2.13) 


(2.14) 
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This result gives an analogous result to the theorem 2.1. Namely, that to find an interval [nL^nL + L] 
such that u{t, x) take the value m on average, one has to choose n such that uo{x) roughly looks like B{m, T) 
for X G [nL, nL + L]. 

In theorem 2.4, the B{m,T) term is a bit mysterious. In the homogeneous case, we characterize the 
location of the level sets in terms of the global-in-time solution to the logistic equation, an ordinary differential 
equation. The following proposition does the same here, making rigorous the ODE-like behavior of u by 
characterizing B{m,T) in terms of solution to a global-in-time analogue to the logistic equation. 

Proposition 2.5. Let the conditions of theorem be satisfied. Let i/jq and /o be the unique positive 
solutions to the eigenvalue problem 


{tpo)xx + fu{x, 0)tjjo = fotjjo, a; G T, 
f t/jgdX = 1. 


Let Tm be the unique time that f (p(Tm)dx = m. Then 

B{m,T) f ifo(x)dx 


lim 


T-s-oo f tp{Tm - T,x)'tpo{x)dx 


= 1 . 


Again we may obtain an easy corollary regarding the movement of the average-level sets Em(t). 

Corollary 2.6. Let the hypotheses of proposition 2.5 be satisfied. For any e > 0, and any m G (0,1) there 
is a time Tq, depending on all parameters, such that if T > Tq, then 

Em(t) C [ug ((^(Tm+e — T)) ,Uq ~ T’))] ■ 

As a part of the proof, we will show that there is a positive constant a > 0 such that, for very large 
negative times, <p{t,x) ~ ae-^°*'tfo{x)- Hence proposition 2.5 and corollary 2.6 imply that there is a constant 
Cm, depending only on m and /, such that 

Emit) ~ Mg ^ (cme“*‘) 


holds for large times t. The constant Cm is, in principle, computable given a specific non-linearity / by 
solving the global in time periodic problem. 

Finally, we prove a result which is analogous to the asymptotic flatness result of [16]. This final result 
shows that u grows exactly as the global in time solution of the periodic problem. 


Theorem 2.7. Let the conditions of theorem 2.f be satisfied. Let tjjQ be as in proposition 2.5. Define 


Then 


Sn = (^J 'tpQix)dx'^ uoinL). 


lim max ||m(T, •) - ipiT^’' + T, ■)\\L°-([nL,nL+L]) = 0. 

T->-oo nGZ 


(2.15) 


Related work 

We will now briefly describe related work. As we have mentioned above, the most closely related work is 
that of Hamel and Roques [16]. In this novel article, the authors present three main results regarding the 
solutions to the homogeneous Fisher-KPP equation (i.e. (1.1) with / depending only on u). First, they show 
that l|Ma;/Ml|L°o(R) —>■ 0 as f tends to infinity. Further, they show that if the initial data decays slower than 
x~°‘F for some ag depending only on the smoothness of /, there are constants Cm and C such that, for t 
large enough. 

Emit) C [mq ^(C'mexp-/'(0)t),Uo ^(Cmexp-/'(0)f)] . 
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We note that a/ = 2 when f(u) = m(1 — u). Lastly, for more general initial data, they obtain the exponential 
order. Precisely, they show that, for any e. 

Emit) C [Mo^(exp-(/'(0)-e)t),Mo^(exp-(/'(0) + e)t)] , 

holds for any t large enough. 

Hamel and Roques obtain their bounds on Emit) by carefully building sub- and super-solutions using the 
function, uoix)e-^ and using the fact that we may bound / below by a second order Taylor approximation 
at zero. Our method involves quite different sub- and super-solutions, and the difference in precision between 
the constant in their bounds and our bounds, in eq. (2.10), can be attributed to the fact that our sub- and 
super-solutions feel the full affect of the non-linearity, /. 

Beyond the work of Hamel and Roques, there has been recent interest in super-linear in time propagation 
phenomena, or accelerating fronts. In a series of articles [7-9], Cabre, Coulon, and Roquejoffre investigated 
the phenomenon of exponential spreading of level sets of solutions to the Fisher-KPP with the fractional 
Laplacian. In both cases, the mechanism is the same: when the data is spread out, ui-,x) evolves according 
to the logistic equation for every x. In this model, the fat tails of the fractional Laplacian are responsible for 
spreading u out in finite time. We note that our techniques can be adapted to yield analogous results when 
uq decays slower than the critical cases investigated by Cabre, Coulon, and Roquejoffre. These accelerating 
fronts have also been observed in biological models with a non-local dispersal operator with fat tails [12,19,23] 
and in a kinetic reaction-transport equation [4]. Finally, they have been conjectured to occur in a reaction- 
diffusion-mutation model with variable motility [3]. 

Plan of the paper 

The results of theorem 2.1 are a special case of theorem 2.4, and, hence, we only prove theorem 2.4. We do so 
in section 3 by creating sub- and super-solutions using the heat equation and eq. (2.12), with carefully chosen 
initial data for both. Then we prove proposition 2.5 in section 4 through an eigenfunction decomposition. 
This characterizes the location of the level sets in terms of the global in time problem eq. (2.11). In section 5, 
we prove theorem 2.7, showing that at large times u locally looks like the global in time solution to eq. (2.11). 
Finally, we show that eq. (2.11) admits a unique global in time solution in section 6. 

Notation 

In order to make the discussion clearer, we adopt the following notation. The constant, C, denotes an 
arbitrary constant independent of time which may change line-by-line. In addition, we occasionally use o(l) 
to mean a constant that tends to zero along with some quantity which we will specify. We will make clear, 
in each usage, the dependencies of C and o(l). 

3 Spreading of the level-average sets 

In this section, we prove theorem 2.4 by creating sub- and super-solutions by decomposing u into a product 
of purely diffusive and purely reactive terms. We utilize the Feynman-Kac formula along with our condition 
on the oscillations, eq. (2.6), in order to control the diffusive term. 

In addition, we state a lemma which will allow us to translate the simple condition on the decay of ug, 
eq. (2.6), to a more useful condition on the oscillations of ug very far to the right. 

Lemma 3.1. Suppose that ug € is strictly decreasing and satisfies 

,. iug)xix) _ 

X^oo Ugix)/ \ogiugix)) 

Let At be some quantity for which there exist eonstants Ci,C 2 > 0 such that At < Then, for any 

choice of cg, we have that 

Ug iuQ^ i\t) ± Cgt) 

lim ---- - = 1. 

i—>-oo At 
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The limit may he taken uniformly for any ci which is bounded above. 

We will use this with the choice A* = B{m,t). To this end, we need the following lemma that guarantees 
the exponential decay of B{m,t) in time, which will be useful in its own right. 

Lemma 3.2. Let B{m,T) be defined as in theorem 2.f. Let fo be the largest eigenvalue of A + fu{x,0). 
Then /o > 0 and there exists positive constants C, independent of m, and Cm, depending on m, such that 

< B{m, T) < Cme-^°^. (3.16) 

Moreover, there exists do > 0 such that if m < 6 q, then 

B{m,T)<Cme-f°^. (3.17) 


We prove these lemmas in section 3.3. 


3.1 An upper bound 

To begin, we define 

Xm{T) = UQ^{B{m,T)), 

and we define a super-solution as U = where satisfies eq. (2.12) and w satisfies 

{ Wt = Wxx + 2 (^7^) Wx, 
uJ(0,x) =max{-2g^,l}. 

Indeed, we notice that z;(0,a;) > uo(x) for every x and that 


Vt - Vxx - fix, v) = -b - tIx'^w - 2ipl’'^Wx - gP^'^^Wxx - fix, 

1 

- - 2ipl^^Wx 


- gF’^Wxx - fix, gF’‘^w) 


Wx 


VTr 


tWx 


= wip^’^' 


fix,gf^’^) fix,Lp^'^w) 


gfT.r. 




> 0 . 


(3.18) 


(3.19) 


The inequality in the last line follows from the fact that w > 1, by the maximum principle, and from the 
fact that fix,s)/s is decreasing in s by eq. (2.3). 

As a result of the above decomposition, we need only show that w does not deviate much from 1 near 
XmiT). To this end, we will use the Feynman-Kac representation of wJ, see e.g. [1]. First, it follows from the 
standard parabolic regularity theory that 

< Co, (3.20) 


where Cq is a constant independent of t, T, and m, see e.g. [20,21]. We define a random process Xt as the 
solution to 



Ir-til] dt + V2dBt, 


(3.21) 


where Bt is a standard Brownian motion. We may then represent w as 


liJ(t,a:) =E"[lF(0,Xt)]. 

Fix ri > 0, to be determined later, and any r 2 . For x > XmiT) — r 2 T, we estimate this at time T as 

wiT, x) = E" [1JJ(0, XT)l{Xr<.-r,T}] + E" [(w(0, Xt) - l)l{;f,>,_,,T}] + E" [l{x,>x-r,t}] . (3.22) 
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First, we may use lemma 3.2 to find a constant C with which we may bound ||?n(0)||oo < . Then, if 

r > 2Co + 2-\/7(), with Co as in eq. (3.20), the first term in eq. (3.22) may be bounded as 


[nJ(0, XT)t{XT<x-rT}\ < Ce*^P {Xt <x-rT}< Ce^“^P - 2CqT < -rTj 


< 


Ce^«^p{V2BT < -2 v^t| 


< 


C 

Vt' 


The last inequality follows from a standard estimate on the tail of a Gaussian. Thus, the first term tends 
to zero as T tends to infinity. The third term in eq. (3.22) can be shown to converge to 1 in the exact same 
manner as above. Notice that these limits did not depend on m. 

Hence, we need only show that the second term converges to zero. To this end we bound the second term 
as 


E" [(nJ(0,XT) - l)t{X,>.-r,T}\ < E" 


< 


'uo{x - riT) - uo{xmiT)y 

Uo{Xm{T)) 

uo{xm{T) - (ri + r 2 )T) - 

B{m,T) 


Using the definition of Xm{T), eq. (3.18), this tends to zero as a consequence of lemma 3.1 and lemma 3.2. 
We notice that this limit can be taken independent of m if m is bounded away from 1 as a consequence of 
the tracking the dependence on m through lemma 3.1 and lemma 3.2. 

Hence we obtain that if a; > Xm{T) — r 2 T, then W{T, x) < 1 + o(l), where o(l) tends to zero as T tends 
to infinity. This yields 

u(r, y) < ^^'-(T, y)wiT, y) < y)[l + o(l)]. (3.23) 

By construction, / <p^’™(r, y)dy = m. Hence, averaging over [a:, a: + L] gives us the upper bound eq. (2.13). 

We note that, to take the limits above, we used only estimates on the Gaussian and lemma 3.1 and 
lemma 3.2. Hence, these limits can be taken independent of m as long as m is uniformly bounded away from 
0 and 1. 


3.2 A lower bound 

The argument here is similar to the argument in section 3.1, above. As before we fix T > 0 and m € (0,1), 
define Xm{T) as in eq. (3.18), and we v{t,x) = ip'^’^{t,x)w{t,x) where (/J^’™(t,x) is defined in eq. (2.12). 
Here, we instead define w as the solution to 


Wf = 


w(0,x) = 


+ 2 ^ 




lin I 


lAo(ai) 



Notice that v(0,a:^) < uo(x) for all x. The computation above, eq. (3.19), shows that 


Vt-Vxx - = Wip^’"" 




‘f 


f{x, Lp^^^w) 

^ipT,m 


where the last inequality follows since n; < 1, as a consequence of the maximum principle, and since /(x, s)/s 
is decreasing for s € [0,1]. 

Again, defining Xt by eq. (3.21), we obtain that 


?i;(t,x) =E"^[w( 0 ,At)]. 

Fix ri > 0, to be determined below, and any r 2 . For any x < Xm{T) + r 2 T, we decompose w{T,x) as 
w{T,x) = E"^ [w(0, AT)l{J(:y>a;+riT}] +E“ [(:U;(0 ,Wt) — l)l{Xt<a:+riT}] pE"^ [l{Xt<a;+r-iT}] ■ 














For X < Xm{T) + 7’2T, we may again choose ri large enough that the first term on the right hand side tends 
to zero as T tends to infinity and that the last term on the right tends to one as T tends to infinity. In 
addition, lemma 3.1 and lemma 3.2 imply that the second term on the right tends to zero as T tends to 
infinity. Hence, we obtain, for any x < Xm{T) + r 2 T, 

u{T,x) > ip'^’'^{T,x)w{T,x) > v?^’™(T,a;)(l-o(l)). (3.24) 

Using that, by construction, j = m, integrating over [x — L,x] yields the lower bound in 

eq. (2.14) for any x < Xm(T) + ?' 2 F. Again, we note that the limit in eq. (3.24) may be taken uniformly for 
m bounded away from 1. 


3.3 Controlling the oscillations 

In this section we prove lemma 3.1 and lemma 3.2. We first prove lemma 3.1, by using an ODE argument to 
compare uq to a sequence of equations which decay at the critical order but with increasingly slower rates. 

Proof of lemma 3.1. We will prove this lemma for Uq ^(At) — cat, though the proof is exactly analogous in 
the case of Uq ^(At) + c^t. Fix e > 0, and define he[x) = e~'^. Notice that 


h' = -- 


2 log(/i) ■ 


(3.25) 


Define 

xq = U(C^(At) and j/o = h~^{Xt). 

We may then shift h to obtain g(x) = h{x — xo + y^). Notice that g[xo) = uo(xo). The condition on the decay 
of At and the fact that uq decays slower than any exponential implies that it()'^(At) — c^t tends to infinity. 
Hence we may choose t large enough, uniformly in ci if ci is bounded above, such that ii x > xq — c^t, then 


0 > u'oix) >-- 


Uo{x) 


\og{uoix)) 


(3.26) 


We note that uq is a supersolution to eq. (3.25) as a consequence of eq. (3.26). Hence uq{x) < g{x) for all 
~ c^t < X < xq. Using this, along with the fact that uo(xo) = g(xo), we obtain 


0 < uo(xo - cst) - uo(xo) < uo(xo - cst) - g(xo) < g(xo - c^t) - g{xo). 


(3.27) 


To conclude, we estimate the last term, g{xQ — c^t) — g{xo). Returning to our function h and manipulating 
the equation, we obtain 

g{xo - cst) - g{xo) = h{yo - c^t) - h{yo) 

= exp{-v^e(2/o - Caf)} - exp{-ye^} 

= exp{-ye^} exp 

It is a simple consequence of Taylor’s theorem that there is a constant C such that 


1 - 


- 1 


1 - a/1 - (cst/yo) < C{c3t/yo). 

Applying this in the equation above, and using that, by construction, = At, we obtain 


g{xo - cit) - g{xo) < At 


exp 


r Cy/ecst I 

I i 


- 1 


Then using that A* < cie we may solve for yo to obtain 

2/0 > \ {c2t - \og{ci)f . 
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Returning to eq. (3.27), and using that uo(xq) = At, by construction, we obtain 


0 < uo(xo - Cat) - At < g(xo - cat) - g(xo) < At 


exp < e 


Ccat 


'c 2 < - log(ci) 

Adding At to to all sides of each inequality and then dividing by At yields 



1 < 


uo{xo - Cat) 

At 


< exp {C't}, 


where C is a universal constant which bounds Ccatjicit — log(ci)). Again, C can be chosen uniformly in 
Cl as long as Ci is bounded above. The result then follows by using the fact that we can choose e as small 
as we would like by taking t larger. □ 


We will now prove lemma 3.2 in order to make lemma 3.1 applicable. The idea of the proof is to use the 
linearized equation to control the growth of B{'m,T). 

Proof of lemma 3.2. First, we verify that /o is positive. To see this, we simply use the Rayleigh quotient 
characterization of /o as 

^ / (/^(x, 0 )V>(x)^ - iVipix)]^) dx 

° ipGH^{T) J,^'tj;{x)'^dx ’ 

where T = [0,L] with the ends identified. One can verify the positivity of /o by testing the equation with 
the function if = 1 and using the positivity of fu{x,0). 

Fix T > 0, and let tpo be the unique eigenfunction associated with fo with Ifyolb = 1- It follows from 
the Krein-Rutman theorem that this function exists and is positive [10]. Notice that Tp{t,x) = a'il)o{x)e^°* 
satisfies 

Pt = + fu{x,Q)Tp> /\lf + f{x,Tp). 

Hence, if a'ifa{x) > B{m,T), we have that Tp(t,x) > (p'^’‘^{t,x) for all x and all t, where (/jT™ jg defined by 
eq. (2.12). Choosing 

B{m,T) 
min^, '00 (a::) ’ 

this implies that 


B{m,T) = ajfyollooe*^ > -f ip{T,x)dx > -f ip'^''^{T,x)dx = m. 

Ta.m.^tpa{x) J J 

Rearranging this inequality yields the lower bound on B{m,T). 

The upper bound is similar, but requires slightly more work. We first obtain a bound for small m by 
using a modulation technique found in [14,15]. Notice that our conditions on / guarantee the existence of 
M such that 

f{x,s)> fu{x,G)s- (3.28) 

holds for all s G [0,1]. Define <f{t,x) = a{t)tpo{x)e^°*. We find conditions on a to make a sub-solution of 
(/?. To this end, we compute 

!£t~ '£.xx~ = ^^ + /o^-a(0o)xxe**-/(x,^) 

= —P + foP - a[/o'0o - fu{x, O)0o]e-^“* - f{x, (p) 

a— — — 

c/ 

= —p + fop-fop + fu{x,0)^-f{x,^) 

< — (p + Mip^'^^ 

a— — 
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Hence is a sub-solution if we choose 


vl + ^ 


^ =-MeW'IWollJ. and a(0) = 


In other words, we may let 


Y{t) = 


B{m,T) 

||'0o||oo 


foS 


-I 1/5 


B{m,TyM[ef°^* - 1 ] + foS 


Hence at time t = T, we may average ip and use that (p sits below to obtain 


B{m,T) 

11'00 11 oo 


foS 


1/5 


efoT^ 


J' p{T, x)dx < J' x)da 


_B{m,TyM[eB^'^ - 1] + fo6_ 

Re-arranging this we see that 

/o<5H(m,T)^e*^^ < m^UoWL +/o<5] . 

If M\\tpo\\^ < foS/2, then we may rearrange the above to obtain 

B{m,Ty <2m^Uo\\Le-^°'^, 


(3.29) 


finishing the claim for small m. 

The argument is somewhat more complicated for general m. We prove it in two steps. First, we show 
that there is Ri, which is uniformly bounded if m is uniformly bounded away from 0 and 1 , such that 
B{m/2,T — Ri) > B{m,T). To see this, we take Cj^rn to be a positive constant, arising from parabolic 
regularity and the smoothness and positivity of /, such that /(x, x)) > Cf^m if the average of p is 

between m/2 and m. Then we choose Ri = ml{2Cf^rn) and we will show by contradiction that this is the 
correct choice of i?i. If not, and instead B{m/2,T — Ri) < B{m,T) holds, then the maximum principle 
gives that 

y = j - Ri,x)dx < j'p'^'"^{T-Ri,x)dx. (3.30) 

Using our choice of Cf^m above, we may integrate the equation for p'^’’^ in space for any t G [T — i?i, T] to 
obtain 

dt j' <^’^’”"(0 x)dx = j' f{x, p'^'"")dx > Cf^m- 

Integrating this in time from T — Ri to T, using the bound in eq. (3.30), and using our choice of i?i yields 
m = ^ p^’'^{T, x)dx > Cf^jnRi + j' p'^’'^{T - Ri,x)dx > y + y, 
which is clearly a contradiction. 

Having proven this fact, we now bootstrap the result for small m to obtain the statement for general m. 
For any m, we may find N such that m/2^ is small enough to apply the work above and obtain eq. (3.29). 
Indeed, we choose N large enough that 


m 


MUoWL < foS/2, 


and as a result obtain that, for any T > 0, 


B 



(3.31) 
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In addition to this, we may iterate the procedure above to find i?i,i? 2 j ■ • ■, Rn such that 


B (^, T - (i?i + • • • + B (^, T - (i?i + • • • + 

holds for any i, and hence, 

B(m, T) < B (^,T - (i?i + • • • Rn)) ■ 

Combining eqs. (3.31) and (3.32) and defining Cm = we obtain 

B(m,T) < - 

where we have defined Cm = Cm2 ~^^-Hn)_ This concludes the proof. 


(3.32) 


□ 


4 Characterization of the speed 


In this section, we will show how to characterize B{m,T) in terms of the the global-in-time solution tp of 
eq. (2.11). In the homogeneous case, there is nothing to prove since (p and ip’^’'^ are equal up to a translation 
in time, by the uniqueness of solutions of ordinary differential equations. In the inhomogeneous setting, 
we need to deal with the fact that (p is not necessarily flat, while (p'^’'^(0,x) = B{m,T). The idea of the 
proof is to run the system for time T/2 and use spectral estimates to show that the ip{Tm — r/2 -|- t) and 
pF’'^{T/2 -I- 1) will be close. From then on, they must remain close. 

Before we begin, we will need one fact about the global in time solution. We delay the proof of this 
lemma until section 6 as it will be crucial in the proof of the well-posedness of eq. ( 2 . 11 ). 


Lemma 4.1. Let tp be any solution of eq. (2.11). 


Then there exist positive constants a and uj such that 


<p{t,x) 

t->-oo ijjQ(^x)ef°* 


and lim 


1 - ipit,x) 

i->-oo 'tj;i[x)e~R'^ 


= UJ 


where ^jJo and fo, defined in proposition 2.5, are the normalized eigenfunction and principle eigenvalue of 
^ + /u( 2 ;, 0 ), and where tpi and fi are the normalized eigenfunction and principle eigenvalue of A — fu{x, 1). 

Now, our first step is to prove the following lemma. 

Lemma 4.2. Suppose that 70 is a C°°(T) function such that || 7 o||oo < Ce~^°^ for some C > 0. 


It = Ixx + /(a;, 7 ), 

7 ( 0 , a;) = 7 o(a;). 

Fix e > 0. Then there is Tq > 0, which can be chosen uniformly if C is bounded above, such that if T > Tq, 
we have that 

i{T/ 2) - tjjo (^J -fotpo^ 

Proof. Define 

r(t, x) = "f{t,x)e^°^'^~*\ 

and we can easily see that F satisfies 


2 


r* - F,, + foT = /(a;,7)e/“(^-‘) = T. 
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Letting ri(t) = f r{t,x)ipo{x)dx, we notice that 

fi = J (/(x, 7 ) - fuix,0)j)i’o{x)dx. 

First, we notice that Fi < 0. In addition, since ipo multiplied by a large constant is a super-solution of F, it 
follows that F is uniformly bounded. We utilize this, along with eq. (3.28), to obtain 


Dividing by Fi, integrating in time, and exponentiating yields 


1 > 


ri(T/ 2 ) 

ri(o) 


^expj-e-^^C"' 


1 -^-^T 

1 — e ^ 


}■ 


Defining 71 (t) = / 'y{t,x)ipQ{x)dx, we may rewrite this inequality as 

|7i(T/ 2) - e*^/^7i(0)| < o(l)e"*'^/^ (4.33) 

Now we let F2(t,a;) = r(t,a;) — Ti{t)tpo{x), and we bound F2. Notice that F2 satisfies the equation 
(r2)t-(F2),,-(A(x,o)-/o)r2 


/(a;,7)e*(^ ‘)-/u(x,0)F - V'o(a;)e*('^ j {f{yn) - fu{y,Q)l{y))'4^o{y)dy 


(4.34) 


Because F 2 is orthogonal to the the eigenspace for eigenvalue /o of A -|- fu{x, 0), then there is a such that 

J [(r 2 ).. + ifuix, 0) - /o)r 2 ] T 2 dx < -a\\T 2 \\l 


Hence, multiplying eq. (4.34) by r 2 and, again, using that F 2 and are perpendicular yields 
^^I|r 2||2 + a:||r 2||2 < [f{x,^) - fu{x,0)-/]T2dx 

< J 7^+'^|r2|fia; = J ri+'^r2da; 

As we noted above, ||r||oo is uniformly bounded in time. Hence, we may write this as 

|||F 2||2 + 2 a||F 2 || 2 <Ce^*(‘-^)||F 2 || 2 , 

where C is a new constant independent of t and T. Solving this differential inequality and using that 
||r 2 ( 0 )|| < C for some constant C, yields constants C > 0 and 0 > 0 such that ||F 2 (T/ 2)||2 < Ce~^^. 

The combination of this inequality with eq. (4.33) yields that 

||7(T/2) - e*^/2^i(0)V’o||2 < C'|7i(T/2) - e*^/2.^i(0)| + e-*^/2||r2||2 < o(l)e-*^/2, 

hnishing the proof. □ 

We now show how to use this lemma to conclude proposition 2.5. Fix m G (0,1) and define the starting 
time S{m,T) as the time such that 

J ip{S{m,T),x)ijjo{x)dx = B{m,T) J ipQ. (4.35) 
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We first show that f (p(S(m, T) + T, x)dx m as T tends to infinity. Then, defining Tm 
such that 



as the unique time 


we will leverage this to show that 


/ ip{Tm - T, x)tlJo{x)dx 
T-s-oo B{m,T) f tpo 


(4.36) 


finishing the proof. 

By parabolic regularity, lemma 4.2, and lemma 3.2, we have that 

MS{m, T) + T/2) - (^^’-(T/ 2)||2 < 


where o(l) tends to zero as T tends to infinity (uniformly for m bounded away from 1). By parabolic 
regularity this may be strengthened to 

MS{m,T) + T/2) - ^^’™(T/2)|U < o(l)e-*^/ 2 _ 

Define r]{t, x) = T) + T/2 + t,x) — ip^’^{T/2 + t, a:) and notice that r] solves 


m - Vxx = c{t, x)ri, 

where c(t,a;) = {f{x,ip^) — f{x,ip'^))/{ip^ — (f'^) < /„(x, 0). Hence, is a super-solution to the 

equation above with initial data larger than r]{0,x). This implies that 

\r]{T/2,x)\ < o(l)a/>(x) < o(l). 


holds for every x. Using the definition of rj, this implies that 

MSim, T) + T, •) - Olloo < o(l). 

By our definition of this implies that 

lim -f ip{S{m,T) +T,x)dx = m, 

T-s-ooJ 

holds, where the limit may be taken uniformly in m if m is bounded away from 1 . 

Now we will show eq. (4.36). Suppose this does not hold. Let 6 > 1, without loss of generality, be such 
that there is a sequence Ti, r 2 , T 3 ,... such that 

f p{Tm-Tn,x)'lpo(x)dx 

lim - - - - —p - > (J. 

n-).oo B{m,Tn) } ipo 


We note that the proof is similar for 0 < 1 and the limit above is less than 9. Using lemma 4.1 and the 
definition of S{m, T), we can choose Nq be such that if n > iVo then 


gT„-T„ 

gS{m,T^) 


and hence, 

T^-T„-5(m,T„) >log(0) >0. 

By parabolic regularity, along with the smoothness and positivity of /, we may find e depending only on m 
and 6 such that f{x, (p{t, x)) > e for all t € [T„ -|- S{m, T„), Tm]- This holds because p is bounded away from 
0 and 1 on this interval. Integrating the equation for p on this time interval yields 

p{Tm) - p{Tn +S{m,Tn)) = [ f {x, p)dx > e\og{6). 

JT„+S(m,T„) 


The average of the first term is m, and the average of the second term converges to m for large n. This 
implies that the left hand side tends to zero, which is a contradiction. This finishes the proof. 
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5 Convergence to the global in time problem 

In this section, we prove theorem 2.7. We will do this in two steps. First we will prove the closeness of u 
and Then we will show that and (p are close when suitably translated in time. To this end, we 

state two lemmas. 

Lemma 5.1. Define as in eq. (2.12), and define Im,T to he the interval [nL,nL + L) containing 

UQ^{B{m,T)), we have that 

lim max ||w(T, •) — •) || , \ = 0. (5.37) 

T-)-oo me(O.l) " " " ^ 

Lemma 5.2. Define S{m,T) as in eq. (4.35). Then, 

lim max \\ip{S{m,T) + t, ■) - ip'^’'^(t,= 0. 

I —^oo t£[T/2,oo\ 


This limit can be taken uniformly for m which is bounded away from 1. 

We will prove eq. (5.37) at the end of this section. The proof of lemma 5.2 follows from our choice of 
S{m,T) and from work in lemma 4.2. Hence, we omit it. We now show how to conclude theorem 2.7 using 
these two lemmas. 

Proof of theorem 2.1. First, denote Impr as [nL,nL + L\. Define T^^ as the unique time such that 


J ip{T^^,x)dx = Sn, 

where Sn is as in eq. (2.15). We show that if T (or equivalently n) is large enough, then T'®" and S{m,T) 
are close, and that we can obtain this in a uniform way for m bounded away from 1. To this end, notice 
that, by the definition of Im,T, 

uo{nL) > B{m, T) > uq^uL + L). 

In addition, by lemma 3.1 and lemma 3.2, we have that 


uo{nL) 

B{m,T) 


1 + o(l). 


(5.38) 


where o(l) tends to zero as T tends to infinity in a uniform way if m is bounded away from 1. Using 
lemma 4.1, we then compute that 


1 + oai = = ( f f _ j f’Qdx _\ 

B{m,T) \ (f ipQ(x)dxf J \f T(S(^^T),x)ipo(x)dxJ 

^ ( /a(^o(x) +o(l))e^°^''"da; \ /_ ! fjpdx _\ ^ 6 -^°'^''" 

Iv {fMx)dxf )\!a[Mx)+o{l))efoSira,T)^,[x)dx) e/oS(-T)^ 

Again, the constant o(l) tends to zero uniformly as long as m is bounded away from 1. This yields the claim 
that T^" — S{m,T) tends to zero as T tends to infinity. 

Fix e > 0. Using lemma 5.1, choose Ti large enough that if T > Ti, then 

max \\u{T, ■) - ^^’™(r, •)||l~( 7 ,„ ,) < e/2. (5.39) 

mG(0,l) 

In addition, choose mi < 1, depending only on / and e, such that if m > mi, then 

^^’^{T,x)>l-e/2, 
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for all a; € T and any T. 

We will consider three cases. First, if n is small enough that > B{mi,T), the maximum 

principle implies that 

^{TS’- +T,x)> (T, a;) > 1 - e/2. 

This, combined with eq. (5.39), implies that 

\\ip(T^”- +T,X)- u{T, ■)\\L°°{[nL,nL+L]) < £, 


for T> Ti. 

The second case is when n is such that j- ,x)dx < B(mi,T) tpo') . This implies that m < mi. 

Indeed, using the dehnition of Sn, we obtain that 


B{mi,T) ■00^ > y (piT^^,x)dx = Sn 


uo{nL) / t/o > B{m,T) 



Thus, B(m,T) < B{mi,T), implying that m < mi. Hence we may find T 2 , depending only on mi, /, and 
e, such that if T > T 2 then T'®" is close enough to S{m, T) that 


||^(yS„ + T, •) - ^{S{m, T) + T, OIU- < e/8. 


In addition, increasing T 2 if necessary, we may apply lemma 5.2 to get that if T > r 2 then 

||^(5(m, T) + T, •) - Olloo < e/16. 

The combination of these two inequalities with eq. (5.39) implies that if T > Ti + T 2 , then 

_ u{T, QnL,nL+L]) < £. 

Finally, we handle the case when neither of these is true. That both of these conditions do not hold, 
along with parabolic regularity, implies that we may find a positive constant 0 < I, which depends only on 
mi and /, such that 

9 J-,x)dx < B{mi,T) . 

Using lemma 3.2 we may find to > 0i depending only on / and mi, such that 

j'ip{T^^,x)dx < B(mi,T-to) ■ 

Following the work above, we then have that B{m, T) < B{mi,T — to). Hence we obtain 

(!-„) = /(I - > /(I - y—Mm)).. > 

The constant C depends only on mi and /, and the last inequality comes from considering a constant 
multiple of as a sub-solution to 1 — on the interval [T — to,T]. This shows that m is bounded 

away from I independent of T. Hence we may argue exactly as in the last paragraph, finishing the proof. □ 

We will now show how to conclude lemma 5.1. The main idea is largely the same: when m is large 
enough, both u and are near 1, and when m is bounded from I, we may take the limit in T uniformly 
in m. 
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Proof of lemma 5.1. Fix e > 0, and we will choose T, independent of m, such that 

To this end, we utilize the work from section 3. 

First, we may use parabolic regularity to find mi < 1 depending only on / and e such that 
1 — < e/2. We may now consider, separately, the cases when m < mi and when m > mi. 

If m > mi, using a sub-solution, v = as in section 3.2, we may choose T, depending only on e 

and /, such that 

w{T,x) > (1 - e)/(l - e/2), 
for any x < Hence, we have that 

1 > u{T,x) > v{T,x) > ip^’'^^{T,x)wt > (1 - = 1 - e. 

for any x < max7^1 ,t- By the maximum principle, we can see that if m > mi. Hence we 

have that, if m > mi 

The case when m < mi follows from our work above. Indeed, since m is bounded away from 1, we may 
choose T such that w and w from section 3.1 and section 3.2, respectively, satisfy 

1 — e/2 < w{T,x), and w{T,x) < 1 + e/2, 

for any x € Im,T- Hence we have that, for any x £ Im,T, 

(l - I) p^^^{T,x) < w{T,x)p^’^{T,x) < u{T,x) < w{T,x)p^’^{T,x) < (l + |) P^’^iT,x), 
which finishes the proof. □ 


6 Well-posedness of the global in time solution 


In this section, we establish the existence and uniqueness (up to translation) of the solution to eq. (2.11). In 
other words, we prove the following proposition. 

Proposition 6.1. There exists a unique global in time solution of eq. (2.11). 

We will prove this in a series of three lemmas. First, we will exhibit the existence of a solution. Afterwards, 
we show that any solution to eq. ( 2 . 11 ) must look grow and decay exponentially as t tends to — oo and oo, 
respectively. Finally, we will use this qualitative property and the sliding method to establish uniqueness. 

Lemma 6.2. There exists a solution of eq. (2.11). 

Proof. We establish existence by considering the limit of family of initial value problems. Define to be 
the solution to the parabolic equation on [T„,oo] x T 

< fp’^(0,x)dx = 1/2, 

I (/j”(T„,x) = 1 /n, 


where T„ < 0 is an unknown. It is easy to check that (nmin'^o) ^'ipo(x)e-^°^* is a super-solution, where 
tpQ and /o are the principle eigenfunction and eigenvalue of A -|- fu{x,0). Thus, we obtain 


1 

2 



||'0o||oo ^-faT„ 
n min i/q 
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This implies that r„ < Clog(n), for some constant C. Hence we may use parabolic regularity, along with 
the compactness of Holder spaces to obtain convergence along a subsequence to a global in time function ip 
which satisfies 

\ipt= ^xx + 

p{0,x)dx = 1 / 2 , 

finishing the proof. □ 

Now we establish that the asymptotics of any solution p to eq. (2.11) look like the solutions to the 
linearized problems around 0 and 1. This amounts to proving lemma 4.1, as stated in section 4. 

Proof of lemma 4.1. We will prove the claim for t tending to — oo. The result for t tending to infinity follows 
by looking at (pit, x) = 1 — p{—t, x) and arguing similarly. Our strategy is to first show that p is bounded 
above and below by a multiple of Then we will appeal to parabolic regularity to obtain the result. 

To this end, we first assume, by contradiction, that e~l°*p(t,x) tends to zero along some sequence 
tending to —oo, for some choice of Xn. By the parabolic Harnack inequality, it must converge to zero uniformly 
in space. Hence, we may assume that, for any e, there exists N such that if n > fV then \\e~ 1°*^ p{tn)\\oo < £■ 
In this case, we may apply lemma 3.2 to find large enough that — tn,x) is a super-solution to 

p on [—tn, 0] X T. Using this inequality at time t = 0, we have that 

\= j X)dx > -j pif), x)dx = 

This is clearly a contradiction, finishing the lower bound. 

Now we will show that there is some C such that e~l°*p{t,x) < C for t < 0. To this end, suppose that 
there is some sequence tending to —oo such that e~l°*^p{tn,Xn) —>■ oo for some sequence a;„. Again, by 
the parabolic Harnack inequality, this implies that ijiine~l°^’^p{tn) —>■ oo. Hence, we may apply lemma 3.2 
and take n large enough that — tn,x) is a sub-solution of p on [—tra,0] x T. Using this inequality 

at time t = 0 yields 

J x)dx < j- p{0, x)dx = i. 

This is a contradiction, giving us the desired upper bound. 

Now, letting = e~l°^p{t,x), we will show that $ converges to atfo. To this end, we may write 

= f ^(t, x)ipo(x)dx and $2 = — '0q4>i. Then as in our work in lemma 4.2, we may see that there are 

positive constants C and 9, independent of time such that 

0 > $1 > -Ce^A‘4.1, 

^||$2|l^ + 20||$2||^<Ce^*‘||<h2||2. 

Solving these differential inequalities and using that $ is bounded for all time gives us that, as t tends to 
—00, 4 >i(t) —>■ a, for some a, and that 4>2(t) —>■ 0. Since $ is bounded away from zero, we get that a > 0. 
This gives us the desired result. □ 

Finally, we use this qualitative information, along with the sliding method, to prove uniqueness of any 
solution to eq. (2.11). 

Lemma 6.3. Any solution to eq. (2.11) is unique. 

Proof. Suppose that p^ and p'^ are both solutions to eq. (2.11) with f p^(0,x)dx = 1/2. Using lemma 4.1, 
we may choose h large enough such that 


p^(t,x) < p'^{t -I- h,x), 
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holds for all t and all x. Then we define 


ho = inf{/i : x) < (p'^{t + h, x)}. 


(6.40) 


Since and (p^ have the same spatial average at t = 0, then hg > 0. We wish to show that hg = 0. To 
this end, let us assume that ho > 0. There are two cases. Using lemma 3.2, we can find ai and a 2 such 
that (p^{t,x) ai'ipo as t tends to —oo that (p^{t,x) a2il’o as t tends to —oo . The first case is that 
ai = e^°^°a2- In other words, ip^ and touch at time t = —oo. Define r]{t,x) = p^(t + ho,x) — p^{t,x) 
and rj satisfies 

_ /(x, p‘^{t + hp, x)) - f{x, p^ (t, x)) 

p'^{t + ho,x) — p^{t + h,x) ^ 

Fix T < 0 to be determined later. Since p'^ ^ and since ai = e-^°^°a 2 , we notice that r]{T,x) < 

o{l)e^°'^ where o(l) tends to zero as T tends to —oo. From this it is easy to see that o(l)'0oe^"* is a super¬ 
solution to equation for r] starting from any time T. This implies that ri{0,x) < o(l). Since this holds 
independent of our choice of T, we get that that 77 ( 0 , x) = 0. 

On the other hand, p'^ is bounded away from 0 and 1 on the time interval [0,/io]- Hence, there exists 
9 > 0 such that f(x, p^') > 0 on the time interval [0, hp]. Integrating the equation for p"^ in time and space, 
along with this fact and the fact that ??( 0 ,x) = 0 , yields 


0 = ^ 77 ( 0 , x)dx = -j' p^{hp, x)dx — -j' 73 ^( 0 , x)dx = -j' p^(hQ,x)d: 


X — m 

ho 


/ r pho p nho 

p'^{ho,x)dx — J p^{0,x)dx = J jf{x,p‘^{s,x))dxds> J 9ds = hp6 > 0. 


The second equality on the second line comes from integrating the equation for p'^ in time and space. The 
above inequality is clearly a contradiction. 

Hence, it must be that, instead, ai < e-^°^°a 2 - Choose e > 0 small enough that a^e'^ < e^°~'^^a 2 - Take 
to to be a very large negative number such that if t < to then 

p^{t,x) < aie‘^ipo{x)e-^°^ 


and 

p‘^{t + ho) > 

Notice then that p'^{t + hp — e,x) > p^{t,x) for all t < tp. In addition, the comparison principle assures us 
that p'^it + ho — t^x) > p^{t,x) for all t > tp as well. This contradicts our choice of hp in eq. (6.40). 

Hence we may conclude that hp must be zero. This implies that p^(t,x) < p'^{t,x) for all t and all x. 
The proof may be applied just as easily to show that p^ < p^ as well. Hence we conclude that p^ = p'^, 
concluding the proof. □ 
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